How to Sum Contributions into the Total 
Charged-Current Neutrino— Nucleon Cross Section 
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The total CC (anti) neutrino- nucleon cross section is usually estimated by the sum of contributions 
from quasi-elastic scattering (QES), single-pion production through baryon resonances (RES), and 
deep inelastic scattering (DIS) with an appropriate scratching the phase space of the RES and DIS 
contributions. However the resulting total cross section is very sensitive to the value of the cut- 
off in invariant mass of the final hadron system produced in RES and DIS. We examine available 
experimental data on the QES and total CC cross sections in order to extract the best-fit value for 
this cut-off. By using the same data set we attempt to adjust the poorly known values of the axial 
mass for QES and RES. 

PACS numbers: 12.15.Ji, 13.15.-|-g, 14.20.Gk, 23.40.Bw, 25.30.Pt 



I. INTRODUCTION 

It is conventional to estimate the inclusive charged and 
neutral current neutrino-nucleon cross sections by the 
sum of contributions from exclusive channels and deep 
inelastic scattering (DIS): 



,tot _ „(Q)ES , 



'vN 



'vN 



(1) 



In the absence of a received model for multi-hadron ex- 
clusive neutrinoproduction, the exclusive contributions 
in Eq. ^ are usually assumed to be saturated by elastic 
(NC case) or quasielastic (CC case) scattering (ES/QES) 
and single-pion production through baryon resonances 
(RES). The exclusive and inclusive (DIS) contributions 
are of the same order of magnitude within the few-GeV 
energy region. Thus, to avoid double counting, the phase 
space of the RES and DIS contributions have to be 
scratched by the conditions W < W^^^ and W > W^\^, 
respectively, where W is the invariant mass of the final 
hadron system in RES or DIS, and W^^^ and W^Jf are 
some parameters. 

The physical basis of this approximation is the concept 
of quark-hadron duality, according to which the result- 
ing total cross section should be essentially independent 
of the specific values of the cutoff parameters if they are 
of the order of the threshold value of W for two-pion pro- 
duction, W27r = Mpf + 2m-^ « 1.2 GeV. In practice, this 
value is too small since the structure functions involved 
into the calculations of the DIS cross section cannot be 
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extrapolated to the two-pion production threshold due 
to the obvious reasons. 



The problem is aggravated by the uncertainties in 
the knowledge of the simplest exclusive contributions: 
the description of the RES reactions is vastly model- 
dependent and, even within a fixed model for RES, both 
RES and (Q)ES cross sections are very sensitive to the 
poorly known shape of the weak axial-vector form fac- 
tors. By adopting the standard dipole parametrization 
for these form factors, their shapes can be described with 
the two phenomenological parameters ("axial masses") 
Af^^^ and M^^^ which, strictly speaking, may be dif- 
ferent and whose experimental values spread within in- 
admissibly wide ranges 0. 



In this study we attempt to fine-tune both the ax- 
ial masses M^^^, Mf^^ and the cutoflts W^J^^, l¥°f 
by fitting available experimental data on the QES with 



AY 



and total CC cross sections for and scat- 



tering off different nuclear targets (converted to the pro- 
ton, neutron, and isoscalar nucleon) as well as the inde- 
pendently measured ratios of these cross sections. For 
the moment, our global likelihood analysis does not in- 
clude the experimental data for the (quasi)elastic reac- 
tions with Ay ^ 0, the reactions of single-meson produc- 
tion and NC induced (exclusive and inclusive) reactions. 
However, the bulk of the data involved into the analy- 
sis is already rather representative and (more important) 
more self-consistent in comparison with the data for the 
single- and multi-hadron neutrinoproduction and the NC 
reactions of any kind. Hence we guess it is sufficient for 
preliminary practical conclusions. 
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II. THEORETICAL MODELS 



B. Resonance single-pion production 



A. Quasielastic scattering 

For the vn p'p and Vp — > p~^n cross sec- 
tions we use the standard result (see, e.g., Ref. 0) 
neglecting the second-class current contributions. For 
the elastic electromagnetic form factors G^" and G^P 
we apply the QCD Vector Meson model by Gari and 
Kriiempelmann |3| extended and fine-tuned by Lomon Q 
to match the current and consistent earlier experimen- 
tal data derived using Rosenbluth separation and polar- 
ization transfer techniques. More explicitly, we explore 
the so-called "GKex(02S)" version of the model advo- 
cated by Lomon. At 4-momentum transfer below 
10 - 15 GeV^ the GKex{02S) model is very close nu- 
merically to the PTD (polarization transfer data based) 
version of the popular "BBA-2003" inverse-polynomial 
parametrization by Budd et al. Q obtained through a 
global fit to the world data on the Sachs form factors, in- 
cluding the results of several more recent measurements. 
Although the up-to-date experiments (see, e.g., numer- 
ous reports in Ref. @) do not contradict to both models, 
we prefer the GKex(02S) model since it meets the re- 
quirements of dispersion relations and QCD asymptotics 
at low and high while the BBA-2003 PTD fit has an 
unphysical behaviour when extrapolated to high (a 
typical drawback of polynomial approximations). 

For the axial and pseudoscalar form factors we use the 
conventional representations 0] 



and 



Fp{Q' 



with Fa(0) = gA = -1.2695 7]. The currently available 
experimental data on the axial mass. Ma = M^^^, show 
very wide spread, from roughly 0.6 to 1.2 GeV/c^ 0. 
Today, it is the main source of uncertainties in the QES 
cross sections. Since the pseudoscalar contribution enters 
into the cross sections multiplied by (me,^,T/.^Ar)^, it is 
substantial for neutrinoproduction of t leptons but small 
for electron and muon production; hence the related un- 
certainty is not important for the present study. 

Since the major part of the experimental data on QES 
obtained for heavy nuclear targets was not corrected for 
nuclear effects, one have to take these into account in 
calculations. We apply the simple Pauli factor since its 
effect for the total cross sections is not essentially dif- 
ferent from that evaluated with the more sophisticated 
approaches. 



In order to describe the single-pion neutrinoproduc- 
tion through baryon resonances we use an extended ver- 
sion of the model by Rein and Sehgal (RS) [1, 0. The 
RS model, being one of the most circumstantial and ap- 
proved phenomenological approaches to calculating the 
RES cross sections, is now incorporated into essentially 
all Monte Carlo neutrino event generators developed for 
both accelerator and astroparticle experiments. Our 
extension Hlf takes into account the final lepton 
mass |l2j and is based upon a covariant form of the 
charged leptonic current with definite lepton helicity. In 
the present calculations, we use the same set of 18th in- 
terfering nucleon resonances with masses below 2 GeV/c^ 
as in Ref. ;8] but with all relevant input parameters up- 
dated according to the current data 0. Significant fac- 
tors (normalization coefficients etc.) estimated in Ref. 
numerically are recalculated by using the new data and 
a more accurate integration algorithm. 

The relativistic quark model of Feynman, Kislinger, 
and Ravndal |l3j adopted in the RS approach unambigu- 
ously determines the structure of the transition ampli- 
tudes involved into the calculation and the only unknown 
structures are the vector and axial-vector transition form 
factors G^'^ (Q^)- -R-^f. Q they are assumed to have 
the form 



G^'-4 (g^) oc 



l/2-r 



K,A 



(2) 



with the "standard" value of the vector mass My — 
0.84 GeV/c^ (that is the same as in the naive dipole 
parametrization of the elastic vector form factor) 
The axial mass Ma = (which was fixed to be 

0.95 GeV/c^ in the basic model) will be a free parame- 
ter in the present study. The integer n in the first ("ad 
hoc" ) factor of Eq. ((SJ is the number of oscillator quanta 
present in the final resonance. 

To compensate for the difference between the SUq pre- 
dicted value (—5/3) and the experimental value for the 
nucleon axial- vector coupling gA, Rein and Sehgal intro- 
duced a renormalization factor Z — 0.75. In order to 
adjust the renormalization to the current world averaged 
value gA = -1.2695 ± 0.0029 (assuming gv = 1) we 
have adopted Z = 0.762. 

Another essential ingredient of the RS approach is the 
nonresonance background (NRB) for which we use the 
ansatz suggested in Ref. [g. The NRB contribution is 
important for description of the existing data on the reac- 
v.in u^jDtt'^, ViiP — ^ u+jDTT^, and 



tions Vfj^n 



p nn ' , ly^n /i pn" , v^p 

also in the RES contribution to the total cross section 
if Wp^f^ < I^c^t^- It is not so obvious for the opposite 
(and by no means unphysical) case, W,,^?^ > , since 

the DIS contribution partially accounts for the NRB. So, 
it would be natural in this case to consider the NRB as 
an additional "free parameter" of the likelihood analysis. 
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However, in this paper we pass over this complication 
and include the NRB contribution into all variants of the 
fit. 

Figure ^ shows the RES contributions into the total 
CC cross sections (divided by neutrino energy) evalu- 
ated by using the extended RS model. In this exam- 
ple, we use M^^^ = 1.08 GeV/c^, the best fit value 
obtained from the recent analysis of the BNL 7-foot 
bubble chamber deuterium experiment (hereafter re- 
ferred to as "BNL-2002" ) based on the total event sam- 
ple of 1.8 M pictures (held two periods of runs in 1976- 
77 and 1979-80). The curves in panels (b) and (c) are 
for the sums of the cross sections for the processes in- 
dicated in the legends. The solid and dashed curves in 
these panels correspond to the calculation with and with- 
out the NRB contributions, respectively. The calcula- 
tions are done for the six values of the cutoff parameter 
= 1.2,1.3,1.4,1.6,1.8,2.0 GeV; clearly the cross 
sections decrease with decreasing the cutoff. 




10 10 

E, (GeV) 



FIG. 1: The total CC single-pion production cross sections 
divided by neutrino energy evaluated with the extended RS 
model using MY'^ = 1.08 GeV/c^ for the six values of the 
cutoff parameter W]^^ = 1.2, 1.3, 1.4, 1.6, 1.8, 2.0 GeV (from 
bottom to top curves in each panel). The solid and dashed 
curves correspond to the cross sections calculated with and 
without the NRB contributions, respectively. 

The next and more substantial drawback of the present 
study is in neglecting the nuclear corrections for the RES 
(as well as for the DIS) contribution. A justification is 



in the fact that these effects were subtracted in a certain 
part of the total cross section data while the necessary 
information is unavailable for another part of the data. 
We intend to remove this drawback in future study. 



C. Deep inelastic scattering 



The DIS CC v^N and V^j^N differential cross sections 
are represented by the standard set of five structure func- 
tions Fi = F^{x,Q^) (see, e.g., Refs. |l^i^): 



2\2 



Y^A, {x,y,E,)F,{x,Q'' 



j2 DIS 

dxdy ~ ^(1 + QV^w)' 

(3) 

where x and y are the usual DIS kinematic variables. The 
coefficient functions Ai are 



Ai 

A2 

A, 



y \ xy + 



= 1-1 



Mnx 



2E„ 



2 J AMnE^ 



(4) 



TO, 



2M^Ei, 

,2 



2MnE^x 



A. = — 



MnE^ 



The functions Fi and F2 are related through the mea- 
surable structure function R = Fl/{2xFi) — ctl/o't, the 
ratio of longitudinal and transverse cross sections in DIS: 



DF2 = 2xFi, D 



1 



1 + R 



1 



(5) 



where v — yE^. In order to satisfy Eq. and, simul- 
taneously, the coUinear parton model (PM) limit that is 
the Callan-Gross relation {F^^ 2xFl^ as ^ 00 
or S) 1), the exact structure functions Fi_2 must be 
related to those in the PM limit, F^^, as 

^^1 = (1 - a + aD) F™, F2 = [a + (1 - a)/D] F^^. 

(6) 

Till the function a = a{x, Q^) is not specified, these rela- 
tions are the most general. There are two simplest limit- 
ing possibilities for a: a = {Fi = F[^, F2 = F^^/D) 
and a = 1 {Fi = T)F[^, F2 = F^^). Our analysis of the 
experimental data described in the next section and test- 
ing many models for the parton density functions (PDF) 
suggests that the "a = 1 model" works quite satisfactory. 
Hereafter we will discuss just this particular case. 

For the structure function R we use a combination of 
two up-to-date parametrizations: inside the nucleon res- 
onance region 1A5 < < 3.9 GeV^ and 0.3 < < 
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5.0 GeV we apply the recent precision result of the Jef- 
ferson Lab Hall C E94-110 Collaboration [11 E^; out- 
side this region we apply the Ra version of the accurate 
6-parameter fit to the world data on R proposed by the 
SLAC E-143 Collaboration 0|. The two parametriza- 
tions are sewn by a 2D B-spline in the boundary of the 
kinematic regions. 

In Fig.|21we show a com par ison of the described model 
with the data from JLab [13 and the results of the ctl 
and (jt separation performed by Dress et al. [23 | and 
based on the data of many older measurements of the ep 
cross sections in the resonance region. The filled bands in 
the figure are obtained by varying within the ranges 
0.18-1 GeV^ (a), 1 - 2 GeV^ (b), 2 - 3 GeV^ (c), and 
3-5 GeV^ (d). 
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FIG. 2: The structure function R = aL/ar vs. obtained 
by the Rosenbluth analysis of the inclusive ep c ross sections 
measured at the JLab Hall C experiment for the 
ranges indicated in the panels. Also shown are the results 
of several earlier experiments on ep scattering in the reso- 
nance region converted to R in Ref. ^ for = 0.8 (a), 
1.1 and 1.4 GeV^ (b). The bands are evaluated by using the 
model for R described in the text and by varying within 
the corresponding ranges; the curves are for the R averaged 
over these ranges. 

Since the JLab fit has been obtained from the data 
on ep scattering, we corrected it to the scattering 
and tested by using the QCD based Altarelli-Martinelli 
equation 0|. In fact, the difference between _R(^'^) and 



jg practically negligible within the relevant kine- 
matic region below the charm production threshold and 
small above the threshold. 

From Eqs. (0J and the exact i^N kinematics it follows 
that 



A4 < 



2MnE^ 



< 



2M 



N 



and 1^5! < 



Mn' 



Due to this suppression and in view of the scale of the 
functions F4 and F5 followed from the NLO pQCD plus 
target mass calculations 16], the A4F4 and Asi^s terms 
in EqJ31) can only be significant near the reaction thresh- 
old [13 ■ Hence the structure functions ^4^5 can be esti- 
mated roughly, by using the approximate relations valid 
in the PM limit with massless quarks: 

2 I 2a; / 2 Id / ' 



and 



Fo F-i 



The PDF contributions into all structure functions are 
divided, in the standard fashion, onto "non charm pro- 
duction" (ncp) and "charm production" (cp) parts: 



^ncp ^ ^ncp^^^ g.) and 



is the Nachtmann 



where xn — 

variable, ^ ~ x n + ni^ / Q'^) is the collinear limit of 
the light-cone variable with massless u, d, and s quarks, 
and rric = 1.3 GeV/c^ is the mass of c quark. The b and 
t quark contributions are neglected. 

We have tested several popular PDF models but in 
this paper we only discuss the results obtained with the 
latest version of CTEQ 6D NLO PDF set with four fla- 
vors (standard DIS scheme, version 6.12, December 14, 
2004) "24^. 

Figure O shows the total DIS f^p, f^rt, T'fj.p, and F^n 
CC cross sections divided by neutrino energy evaluated 
with the CTEQ 6D NLO PDFs for the five values of the 
cutoff parameter W^^f = 1.2, 1.4, 1.6, 1.8, and 2.0 GeV; 
Clearly, the cross sections increase with decreasing of the 
cutoff. Since the CTEQ 6D PDFs cannot be extrapolated 
to the exact kinematic boundaries, we have to freeze 
below some value Q ^ . In Fig. |31 this value varies within 

the range 0.6 to 1.0 GeV^ and the widths of the bands 
reflect the corresponding variations of the DIS cross sec- 
tions. The Qf dependence is in general nonmonotonic 
and diminishes with increasing the cutoff value. In the 
present likelihood analysis, the variable is freezing be- 
low ^ 0.8 GeV^. The error introduced by this approx- 
imation is estimated to be less than 1-2% that is small 
in comparison with the uncertainties of the experimen- 
tal data and indetermination in other phenomenological 
parameters. 
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FIG. 3: The total DIS ffj.p, Vf^n, Vfj.p, and Vfj,n CC cross sec- 
tions divided by neutrino energy evaluated with the CTEQ 6D 
NLO PDFs for the five values of the cutoff parameter W^f — 
1.2 to 2.0 GeV from top to bottom with steps of 0.2 GeV. The 
widths of the bands correspond to variation of the freezing pa- 
rameter Q? from 0.6 to 1.0 GeV^ 



III. DATA SET 

We have examined and classified all available exper- 
imental data on the QES and total CC i^N and VN 
cross sections as well as independently measured rela- 
tive quantities like the ratios <Jvn/<^up, '^vnl'^vpi '^vpl'^vpi 
and so on. Published results from the relevant experi- 
ments at ANL ilJIiJglJlL BNlJgl^ 



FNAL "34 
and IHEP 





CERN 



55, 56 

p 

loa 194 1951 1961 1971 1981 iQiJIJ are 

included dating from the end of sixties to the present day, 
covering v^, 17^, i/e, and v^, beams on a variety of hydro- 
gen and nuclear targets, with energies from the thresh- 
olds to about 350 GeV. A detailed description of our 
database will be published elsewhere. Here we briefly 
depict the most important points. 

Not all the collected data are involved into the analysis. 
We excluded from the fit: 

• the experimental results which are undoubtedly ob- 



solete, superseded or reconsidered (due to increased 
statistics, revised normalization, etc.) in the pos- 
terior reports of the same Collaborations; 

• the datasets which are a transformation of the oth- 
ers derived from the same experimental samples 
(for instance, we used either the cross section u 
or the "slope" ajEy measured in the same experi- 
ment); 

• the cross sections, slopes, and ratios averaged over a 
wide energy range when the energy-binned dataset 
is available. 

We quenched a wish to reject the results seeming self- 
contradictory or being in obvious disagreement with the 
major dataset. A few exceptions and particular cases will 
be expounded in the next section. 

If only the bounds of the energy bin were available, 
we either averaged the data (and the relevant calculated 
quantity) over the bin or estimated the mean energy from 
the (anti)neutrino beam spectrum (when the necessary 
information was accessible from the original paper or an- 
other description of the experiment). The statistical and 
systematic errors of the data were always summed up 
quadratically. 



IV. LIKELIHOOD ANALYSIS 

The four above-mentioned parameters Af^^^, Af^^^, 
W^^f ^, and Wc°t^ involved into the merging of the QES, 
RES, and DIS contributions (Sect. ^ were fitted to 
the described data by using the CERN function mini- 
mizati on a nd error analysis package "MINUIT" (version 
94.1) 100]. In order to test validity of the dataset and 
the fitting procedure, we have examined many variants of 
1-, 2-, 3-, and 4-parameter fits taking care of getting the 
correct correlation coefficients printed by MINUIT. Illus- 
trative results of this analysis are listed in TablesUl Inland 
^3 together with the obtained values of per number 
of degrees of freedom (NDF). The number of significant 
digits shown in the last columns of Tables and IIIII is 
more than needed; we only keep these to clarify the 
minima. 

The first column in each table is for designation of 
different exercises of the fit. The numbers in bold- face 
correspond to the fixed trial values of the parameters 
used as inputs. The errors of the output parameters cor- 
respond to the usu al on e-standard-deviation (Icr) errors 
(MINUIT default) (10l| . 

Visualization of the results is shown in Figs.^l j^for the 
B3 variant of the fit which is preferable in our opinion. 

Fi£ureJlshows the QES data from Refs. 25^ 26, 27, 3j, 
SlIMlsfMSH^IziElillMilllEzllQi to- 
gether with the B3 best fit to the full set of the data satis- 
fying the criteria described in Sect. IIIII fNDF = 670^3 — 
667). The FNAL 1984 data points from Ref. [FJl (neon- 
hydrogen target) are shown here only for a comparison. 
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TABLE III: Three- and four-parameter fits witfi tfie corre- 
sponding per NDF (= 667 and 666, respectively). Trial 
parameters are bold-faced. 



Fit 


Mj 


Ma 


^RES 

cut 


* ^ cut 




NDF 




(GeV) 


(GeV) 


(GeV) 


(GeV) 




0.93 ±0.01 


1.06 ±0.07 


1.27 ±0.01 


1.35 


1.542 




0.95 ±0.01 


1.09 ±0.08 


1.28 ±0.01 


1.40 


1.481 


A3 


0.98 ±0.01 


1.13 ±0.08 


1.30 ±0.02 


1.45 


1.461 




0.98 ±0.02 


1.00 ±0.03 


1.55 ±0.04 


1.50 


1.467 




0.98 ±0.02 


1.05 ±0.03 


1.69 ±0.07 


1.55 


1.463 


B3 


0.98 ±0.02 


1.02 ±0.04 


1.50 ±0.02 


1.468 




0.98 ±0.01 


1.01 ±0.03 


1.50 


1.50 ±0.00 


1.468 




0.98 ±0.02 


1.02 ±0.04 


1.55 


1.51 ±0.02 


1.466 




0.98 ±0.02 


1.03 ±0.04 


1.60 


1.53 ±0.02 


1.464 




0.98 ±0.02 


1.04 ±0.04 


1.65 


1.54 ±0.02 


1.463 




0.98 ±0.02 


1.05 ±0.04 


1.70 


1.55 ±0.02 


1.463 


C3 


0.98 ± 0.02 


1.06 ±0.04 


1.75 


1.56 ±0.02 


1.464 




0.98 ± 0.02 


1.06 ±0.04 


1.80 


1.57 ±0.02 


1.464 




0.98 ± 0.02 


1.06 ±0.04 


1.85 


1.58 ±0.02 


1.465 




0.98 ± 0.02 


1.07 ± 0.04 


1.90 


1.59 ±0.02 


1.467 




0.98 ±0.02 


1.07 ±0.04 


1.95 


1.59 ±0.02 


1.468 




0.98 ±0.02 


1.07 ±0.04 


2.00 


1.60 ±0.02 


1.469 


D3 


0.98 ±0.01 


1.08 


1.73 ±0.26 


1.57 ±0.05 


1.464 




0.8 


1.10 ±0.08 


1.28 ±0.01 


1.38 ±0.01 


1.631 




0.9 


1.11 ±0.08 


1.29 ±0.02 


1.42 ±0.01 


1.497 


E3 


1.0 


1.13 ±0.08 


1.30 ±0.02 


1.46 ±0.01 


1.463 




1.1 


1.17 ±0.06 


1.31 ±0.02 


1.50 ±0.01 


1.551 




1.2 


1.19 ±0.06 


1.32 ±0.02 


1.53 ±0.01 


1.783 


A4 


0.98 ±0.02 


1.13 ±0.08 


1.30 ±0.02 


1.45 ±0.01 


1.463 



value obtained from the global B3 fit. The grey bands 
show the standard deviation from the best-fit cross sec- 
tions due to the error of 0.02 GeV/c^ in determination of 
^jQES^ Note that the best-fit value of M'^'^ is in agree- 
ment with that obtained by a single-parameter fit to the 
QES data only, uf'^ = 0.94 ± 0.04 GeV/c^. 

The obtained value of M^^^, being lower, does not 
contradict to the latest (still preliminary) result by the 
K2K experiment [lO^ 



1.06 ± 0.03 (stat.) ± 0.14 (syst.) GeV/c^ 



It is however essentially below the value of 1.1 GeV/c^ 
used in the recent atmospheric neutrino osci llatio n anal- 
ysis of the Super-Kamiokande I experiment |l04j . 

In Figs. Eland El we collect the main subset of the ex- 
perimental data on the total CC cross sections and their 
slopes for an isoscalar nucleon (hereafter denoted by AQ 
from Refs^_^28J3^08J,l^^ 72, 73, 9^ 

and IMIsalHlMMMlHll^^ 



TABLE I: One-parameter fits with the corresponding per 
NDF = 669. Trial parameters are bold-faced. 



Fit 






H/<5f 




A 


IN UX' 




(GeV) 


(GeV) 


(GeV) 


(GeV) 




0.8 


1.08 


1.47 ±0.02 


1.67 




0.9 


1.08 


1.50 ±0.01 


1.52 


Al 


1.0 


1.08 


1.53 ±0.02 


1.47 




1.1 


1.08 


1.56 ±0.01 


1.56 




1.2 


1.08 


1.58 ±0.02 


1.80 




0.8 


2.0 


1.47 ±0.01 


1.68 




0.9 


2.0 


1.52 ±0.01 


1.52 


Bl 


1.0 


2.0 


1.58 ±0.01 


1.47 




1.1 


2.0 


1.64 ±0.01 


1.56 




1.2 


2.0 


1.71 ±0.01 


1.82 



TABLE II: Two-parameter fits with the corresponding per 
NDF = 668. Trial parameters are bold-faced. 



Fit 


^QES 


^RES 








NDF 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


A2 


0.99 ±0.02 


1.08 


1.53 ±0.02 


1.469 


B2 


0.96 ±0.01 


1.31 ±0.01 


1.4 


1.484 




0.8 


1.11 ±0.08 


1.28 ±0.01 


1.4 


1.631 




0.9 


1.10 ±0.08 


1.28 ±0.01 


1.4 


1.499 


C2 


1.0 


1.09 ±0.07 


1.28 ±0.01 


1.4 


1.493 




1.1 


1.08 ±0.07 


1.28 ±0.01 


1.4 


1.635 




1.2 


1.07 ±0.06 


1.29 ±0.01 


1.4 


1.943 




0.8 


0.91 ±0.04 


1.41 ± 0.02 


1.649 




0.9 


0.97 ±0.04 


1.46 ±0.02 


1.507 


D2 


1.0 


1.03 ±0.04 


1.51 ±0.02 


1.468 




1.1 


1.07 ±0.04 


1.55 ±0.02 


1.558 




1.2 


1.10 ±0.04 


1.59 ±0.02 


1.801 




0.8 


1.00 ± 0.04 


2.0 


1.53 ±0.02 


1.651 




0.9 


1.04 ± 0.04 


2.0 


1.57 ±0.02 


1.505 


E2 


1.0 


1.08 ± 0.04 


2.0 


1.61 ±0.02 


1.469 




1.1 


1.11 ±0.04 


2.0 


1.65 ±0.02 


1.566 




1.2 


1.14 ±0.04 


2.0 


1.68 ±0.02 


1.814 




0.8 


1.08 


1.28 ±0.01 


1.38 ±0.01 


1.629 




0.9 


1.08 


1.77 ±0.09 


1.56 ±0.02 


1.500 


F2 


1.0 


1.08 


1.70 ±0.14 


1.57 ±0.03 


1.465 




1.1 


1.08 


1.34 ±0.03 


1.49 ±0.01 


1.551 




1.2 


1.08 


1.37 ±0.04 


1.53 ±0.02 


1.787 



They are not included into the fit since were obtained 
by a recalculation from the DIS data (included into the 
fit, see Fig. O by using a prescription given in Ref. 'si'l 
and the errors for these points were estimated approxi- 
mately. In order to facilitate comparison, the data points 
for the experiments performed with the nuclear targets 
different from D^and Ne-H2 are converted to a free 
nucleon ta rget 1 102ll. The nuclear effects for the deu- 
terium [2^ I23, 13ll l47l Is^ , neon- hydrogen [sJl and aver- 
aged iron data [53 (shown in Fig. ^ by filled rectangles) 
were subtracted by the authors of the experiments. The 
curves are calculated with A'ff'^ = 0.98 GeV/c^, the 



S IMillBlHlzIlzalmH il Siipi, 

respectively. The majority of these data is included into 
the fit. The curves and bands show the QES, RES, and 
DIS contributions (Fig. O and their sums (both figures) 
calculated with the best-fitted values of uf^'^, Mj^^^, 
and Wp^^^ = Wg^t^ (the latter equality is the restriction 
used in the B3 fit). 
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* Mann et al., ANL 1973 (D^) 
o Barishetal.,ANL 1977(DJ 
o Baker et al, BNL 1981 (D^) 

o Kitagakietal., FNAL1983 (D^) 

* Asratyan et al, FNAL 1984 (Ne-H^) 
Suwonjandee, NuTeV 2004 (Fe) 

* Auerbach et al, LSND 2002 fC") 
•1 Young, CERN 1967 (CF fir) 

Budagov et al, CERN 1969 (C,H,) 



E, (GeV) 



■ Bonetti et al. GGM 1977 (CFfir) 

4 Armenise et al, GGM 1979 {Cflg-CFflr) 

T Pohl et al, GGM 1979 {Cfl^-CFfir) 

* Alla.iia et al, BEBC 1990 (D^) 

a Makeev et al, IHEP SKAT 1981 (CFflr) 

^ Belikovetal, IHEP 1982 (Al) 

□ Belikov et al, IHEP-ITEP 1985 (Al) 

a Grabosch et al, IHEP SKAT 1988 (CFfir) 

a Brunner et al, IHEP SKAT 1990 (CF fir) 

B Ainmosov et al, IHEP SKAT 1992 (CFfir) 



S W 



Barish et al, ANL 1979 
Baker etal., BNL 19S2 
Benvenuti et al, HPWF 1974 
Barish etal, CF 1975 
Barish etal, CFR 1977 
Asratyan et al, FNAL 1984 
Budagov et al., HLBC (F) 1969 
Budagov et al, HLBC (P) 1969 
EichtenetaL, GGM 1973 
Erriquez et al, GGM 1979 
Ciampolillo et al, GGM 1979 
Baranov et al, IHEP SKAT 1979 




* Blietschau et ai, GGM 1978 [a'-' ^ ] 



(GeV) 



NuTeV 2004 
CERN 1969 159 
GGM 1979 




FIG. 4: Total QES cross sections measured by the exper- 
iments ANL 1969 ANL 1973 ANL 1977 \ 
BNL 1981 El, FNAL 1983 El> FNAL 1984 
LSND 2002 :57|, CERN 1967 
, GGM (Gargamelle) 1977 
BEBC 1990 871, IHEP SKAT 1981 
IHEP 1982 93], IHEP-ITEP 1985 94], IHEP SKAT 1988 
IHEP SKAT 1990 and IHEP SKAT 1992 Both sta- 
tistical and total errors are shown for the earliest low-energy 
data of CERN IH^ (excluded from the fit) and for the most 
current high-energy data of NuTeV IssI . The filled rectangles 
are for the NuTeV data (with the total error) averaged over 
the wide energy range 30 to 300 GeV. The data for nuclear 
targets (indicated in the parentheses in the legend) are 
converted to a free nucleon. The curves are for the QES cross 
sections calculated with the value of M^^^ = 0.98 GeV/c^ 
obtained from the global B3 fit (see text). The narrow grey 
bands show the standard la deviation from the best-fit 



FIG. 5: Total CC cross sections for z/,, and i/f^ scattering off an 
isoscalar nucleon measured by the experiments ANL 1979 :23|, 
BNL 1982 0JIPWF 1974 CF 1975 "sSj, CFR 1977 0), 
FNAL 1984~0|, HLBC 1969 Tfreon, 1963-64 and propane, 
1967 runs) jej GGM 1973 IgZ], GGM 1979 III, and 
IHEP SKAT 1979 Also shown are the v^N and V^N 

cross sections measured by the GGM 1978 experiment 
The antineutrino data are scaled with a factor of 0.1 for better 
visualization. The curves and bands show the cross sections 
calculated with the best-fitted values of M^"^^, M™^, and 



(see text and legend in Fig. |SJ . 



Figure [71 accumulates the data on the cross section ra- 
tios (Ju^N / (^u..N and (Ju,^nJoun (isoscalar target) accord- 

m a 154 ^ ^ m £5., m iii 

il Ea El M iil- The paper (eg su- 
persedes the earlier publications of the Gargamelle Col- 
laboration |£^] (shown by filled rectangle) and |6^. The 
major part of these data is obtained from the cross sec- 
tions measured in the same experiments. The other, like 
the recent NuTeV result (HH, correspond to a wide en- 
ergy range with no indication of the mean energy. Due to 
these and similar reasons all these data are excluded from 



the analysis and only shown here for a comparison with 
the result of the global B3 fit. The cross section ratios 

Jp El El US m H 

96] are shown in Fig. |H1 
The results of Refs. ^33j 4Q, , Si, .22., are excluded from 
the fit due to the reasons mentioned in Sect. IIIII The 
near-threshold point from Ref. E3 is removed since its 
deviation from the theoretical prediction is unphysically 
high. Six panels of Fig.|H|show the data of different kinds 
from Refs. [2lElE3|,|42ki2aElElE3- Almost aU data 
points participate in the analysis. The curves and bands 
in Figs. [3 IHl and |51 are calculated with the parameters 
obtained from the global B3 fit (see legend in Fig.EJ. 
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QES:M\ =0.9S±0.02GeV 
RES: Mf^ = 1.02 + 0.04 GeV 
= 1.50 ±0.02 GeV 
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(GeV) 



FIG. 6: Slopes of the total CC cross sections for and 
77^ scattering off an isoscalar nucleon measured by the ex- 



periments ANL 1979 l2g, BNL 



CFRR 1968 12: 
CFR 1977 U 
FNAL 1983 l4i 
CCFRR 1984 \ 
NuTeV 2005 (5; 
BEBC 1979 F 
GGM 1981 IZ! 
BEBC 1984 
ITER 1978 





BNL 1982 
CF 1975 



HPWF 1974 
CFRR 1981 Es 
HBF 1983 m 
, CCFR 1990 |53 
HLBC 1969 JeC 
GGM 1979 IH p 
BEBC 1982 E 
CDHS 1987 85], CHARM 1988 [S^, IHEP- 
IHEP SKAT 1979 Hi, IHEP-ITEP 1979 il, 



FNAL 
FNAL 
CCFR 
GGM 
CHARM 1981 
CDHS 1983 



1982 
1984 
1997 
1973 



cmVCeV (for ly^N) and (0.334 ± 
GeV (for v^N) obtained by the Parti- 



and IHEP-JINR 1996 [991 . The data points with horizon- 
tal error bars are for the slopes averaged over the wide en- 
ergy ranges; these do not participate in the fit and the cor- 
responding energy binned data (included into the fit) are 
shown in Fig. |K| The curves and bands show the QES, 
RES, and DIS contributions and their sums calculated with 
the best-fitted values of the parameters depicted in the leg- 
end in top panel. The averaged values over all energies 
(0.677 ± 0.014) X ^ 
0.008) X 10"^* cm 

cle Data Gr oup llORll from the data by the experiments in 
Refs. |5a . l53ll53 . l85| are also shown for a comparison (straight 
lines). 



V. SUMMARY 

Our analysis of the world neutrino data on the QES 
and total CC cross sections yields several thought- 
provoking conclusions. As is seen from Tables HI ITU and 



FIG. 7: The ratio (j-^n/o'v^n for an isoscalar nucleon mea- 

HPWF 1974 
HPWF 1976 I39|, CFR 1977 
NuTeV 2005 56], GGM 1973 
BEBC 1979 IzJ, CHARM 1981 
(revised according to Ref. |82 
CDHS 1987 m, CHARM 1988 
, IHEP-ITEP 1979 [H, and IHEP- 
The ratio (Tu^jv/ct^^jv reported in the three 



sured by the experiments HPWF 1973J3; 
WHPC 1974 
CCFR 1997 
CDHSB 1977 
BEBC 1982 
BEBC 1984 
IHEP-ITEP 
JINR 1996 |9 

publications of the Gargamelle collaboration |63l l66l |68| is 
also shown. The curve and band are calculated with the 
same values of the fitted parameters as in Fig. |S| 




mil in all variants of the fit there is a distinct minimum of 
for M^^^ around the "canonical" value of 1 GcV/c^ 
with deviations < 2%. This is mainly an effect of the 
QES data subset whose exclusion from the analysis would 
lead to an essential increase of Af^^^ for all variants (for 
example, in the B3 and A4 fits M^^^ becomes equal to 
1.13 ± 0.03 and 1.17 ± 0.03 GeV/c^, respectively). 

The situation with the AI^^^ best-fit value is less def- 
inite: in different variants of the fit it fluctuates from 
about 1.00 to about 1.15 GeV/c^. This spread com- 
prises the BNL-2002 resuhs for M^^^ TB'l obtained with 
different approaches and does not contradict to the ex- 



act equality M 



RES 
A 



M 



QES 



However, the 3- and 4- 



parameter fits favour the case M^^^ > M^^^. Our "fa- 
vorable" B3 variant of the fit yields the following values 



9 



•2 10 

0^ 



u Barish et al, ANL 1979 
® Jacques et al, BNL 1981 
» Baker etai, BNL 1982 
" Hanlon et al, FNAL 1980 
o Myatt & Perkins, HLBC 1971 
H Lerche et al, GGM 1978 
H Armenise et al, BEBC 1981 
m Parker et al, BEBC 1984 
o Allasia et al, BEBC 1984 
H Kayis-Topaksu et al, CHORUS 2003 
a Brunner et al, SKAT 1989 




n Sciulli, FNAL-Michigan 1978 
B Efremenko et al, FNAL 1979 
' Hanlon et al, FNAL 1980 
* Asratyan et al, FNAL 1984 
o Erriquez et al, GGM 1979 
m Abramowicz et al, CDHS 1984 
■ Parker et al, BEBC 1984 
o Allasia et al, BEBC 1984 
B Brunner etal, SKAT 1989 





1.8 




1.6 




1 4 




1.2 




J 








8 


til* 


6 




0.4 


~o 


0.2 




0^ 




0.8 




0.6 




0.4 




0.2 




Q 


=^ 






2.S 




2.6 




2.4 




2.2 




2 




1.8 




1.6 




1.4 




1.2 




J 




0.8 




0.6 




0.4 




0.2 









10 




^ Fartourakis et al, BNL 1980 

• Bahtr et al, BNL 1982 

• Abramowicz etal, CDHS 1984 
B Parker a al, BEBC 1984 

Allasia etal., BEBC 1984 

• Aderhalz el at . BEBC 1 9S6 




i Barish el al, FNAL 1978 
m Parker et al, BEBC 1984 
- Aderholz el at. BEBC 1 9S6 




m Barish el al, ANL 1979 
( QES events removed) 




(d) 
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Abramowicz el al, CDHS 1984 
Parker et al, BEBC 1984 
Aderholz el at. BEBC 1986 




(GeV) 
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FIG. 8: The ratios C7vn/crvp and (TT7n/o"77p measured 
by the experiments ANL 1979 El BNL 1981 [3||, 
BNL 1982 Is^^NAL-Michigan 1978 j40|. FNAL 1979 E|, 
FNAL 1980 H, FNAL 1984 H, HLBC 1971 pi, 
GGM 1978 Jp, GGM 1979 1^, BEBC 1981 M, 
CDHS 1984 137^ BEBC 1984 iSaiSil, CHORUS 2003 H, 
and IHEP SKAT 1989 H. The data point of CDHS 1984 is 
recalculated in Ref. [SOjl from the ratio cri7p/crj7jv. The curves 
and bands are calculated with the same values of the fitted 
parameters as in Fig. |S| 



FIG. 9: (a), (b), (c), (d) - the slopes of the u^p, u^n, 
i7^p, and I7^n total CC cross sections measured by the ex- 
periments BNL 1980 [13, BNL 1982 CDHS 1984 
BEBC 1984 mil, and BEBC 1986 [831; (e) - the ratios 
(7vn/o'up (with quasielastic events removed) and crvp/crvp mea- 
sured by the experiments ANL 1979 I2I FNAL 1978 jH, 
BEBC 1984 and BEBC 1986 |8^ : (f) - the ra- 

tios Ovp/a-^N and a-pja-gN measured by the experiments 
CDHS 1984 73], BEBC 1984 [si, and BEBC 1986 0. The 
curves and bands in all six panels are calculated with the same 
values of the fitted parameters as in Fig. |H] 



of the axial masses: 

^^QES ^ Q gg _^ Q Q2 GeV/c^ 

and 

= 1.02 ± 0.04 GeV/c^ 

The shape of the total and (all the more so) differential 
vN and vN cross sections is very sensitive to the values of 
the cutoff parameters W^^f^ and W^^^^t^. From our anal- 
ysis we have to conclude that these parameters cannot 
be fine-tuned with the confidence level sufficient for the 
current and future experiments for neutrino oscillations 
and related phenomena. However, the most worth-while 
versions of the fit indicate that W^^^^ « W^™? must be 
essentially above the value of 1.4 GeV approved in the 
data processing of many accelerator and astrophysical 
neutrino experiments. The outcome of the B3 fit is 

Wf^^ = = 1.50 ± 0.02 GeV. 



Being considered deliberately, such a high value of the 
cutoff parameter for DIS puts forward the difficult prob- 
lem of a correct accounting for the reactions of exclusive 
multi-hadron neutrinoproduction and coherent neutrino- 
nucleus scattering. 

Finally we have to note that the above conclusions are 
only valid for the theoretical models of the RES reactions, 
DIS structure functions and PDF, as well as the approxi- 
mations and simplifications (sometimes risky) adopted in 
the present analysis. Investigation of alternative models, 
a more accurate treatment of the nuclear effects, and in- 
corporation of additional experimental data is the matter 
of a forthcoming work. 
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